QTQTN motif upstream of the furin-cleavage site plays a key role in SARS-CoV-2 infection and pathogenesis.

The furin cleavage site (FCS), an unusual feature in the SARS-CoV-2 spike protein, has been spotlighted as a factor key to facilitating infection and pathogenesis by increasing spike processing. Similarly, the QTQTN motif directly upstream of the FCS is also an unusual feature for group 2B coronaviruses (CoVs). The QTQTN deletion has consistently been observed in in vitro cultured virus stocks and some clinical isolates. To determine whether the QTQTN motif is critical to SARS-CoV-2 replication and pathogenesis, we generated a mutant deleting the QTQTN motif (ΔQTQTN). Here, we report that the QTQTN deletion attenuates viral replication in respiratory cells in vitro and attenuates disease in vivo. The deletion results in a shortened, more rigid peptide loop that contains the FCS and is less accessible to host proteases, such as TMPRSS2. Thus, the deletion reduced the efficiency of spike processing and attenuates SARS-CoV-2 infection. Importantly, the QTQTN motif also contains residues that are glycosylated, and disruption of its glycosylation also attenuates virus replication in a TMPRSS2-dependent manner. Together, our results reveal that three aspects of the S1/S2 cleavage site-the FCS, loop length, and glycosylation-are required for efficient SARS-CoV-2 replication and pathogenesis.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants.

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for next-generation sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called 'Tiled-ClickSeq', which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5'UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants.

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for Next-Generation Sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called 'Tiled-ClickSeq', which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5'UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

Loss of furin cleavage site attenuates SARS-CoV-2 pathogenesis.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-a new coronavirus that has led to a worldwide pandemic1-has a furin cleavage site (PRRAR) in its spike protein that is absent in other group-2B coronaviruses2. To explore whether the furin cleavage site contributes to infection and pathogenesis in this virus, we generated a mutant SARS-CoV-2 that lacks the furin cleavage site (ΔPRRA). Here we report that replicates of ΔPRRA SARS-CoV-2 had faster kinetics, improved fitness in Vero E6 cells and reduced spike protein processing, as compared to parental SARS-CoV-2. However, the ΔPRRA mutant had reduced replication in a human respiratory cell line and was attenuated in both hamster and K18-hACE2 transgenic mouse models of SARS-CoV-2 pathogenesis. Despite reduced disease, the ΔPRRA mutant conferred protection against rechallenge with the parental SARS-CoV-2. Importantly, the neutralization values of sera from patients with coronavirus disease 2019 (COVID-19) and monoclonal antibodies against the receptor-binding domain of SARS-CoV-2 were lower against the ΔPRRA mutant than against parental SARS-CoV-2, probably owing to an increased ratio of particles to plaque-forming units in infections with the former. Together, our results demonstrate a critical role for the furin cleavage site in infection with SARS-CoV-2 and highlight the importance of this site for evaluating the neutralization activities of antibodies.

Furin Cleavage Site Is Key to SARS-CoV-2 Pathogenesis.

SARS-CoV-2 has resulted in a global pandemic and shutdown economies around the world. Sequence analysis indicates that the novel coronavirus (CoV) has an insertion of a furin cleavage site (PRRAR) in its spike protein. Absent in other group 2B CoVs, the insertion may be a key factor in the replication and virulence of SARS-CoV-2. To explore this question, we generated a SARS-CoV-2 mutant lacking the furin cleavage site (ΔPRRA) in the spike protein. This mutant virus replicated with faster kinetics and improved fitness in Vero E6 cells. The mutant virus also had reduced spike protein processing as compared to wild-type SARS-CoV-2. In contrast, the ΔPRRA had reduced replication in Calu3 cells, a human respiratory cell line, and had attenuated disease in a hamster pathogenesis model. Despite the reduced disease, the ΔPRRA mutant offered robust protection from SARS-CoV-2 rechallenge. Importantly, plaque reduction neutralization tests (PRNT 50 ) with COVID-19 patient sera and monoclonal antibodies against the receptor-binding domain found a shift, with the mutant virus resulting in consistently reduced PRNT 50 titers. Together, these results demonstrate a critical role for the furin cleavage site insertion in SARS-CoV-2 replication and pathogenesis. In addition, these findings illustrate the importance of this insertion in evaluating neutralization and other downstream SARS-CoV-2 assays. IMPORTANCE: As COVID-19 has impacted the world, understanding how SARS-CoV-2 replicates and causes virulence offers potential pathways to disrupt its disease. By removing the furin cleavage site, we demonstrate the importance of this insertion to SARS-CoV-2 replication and pathogenesis. In addition, the findings with Vero cells indicate the likelihood of cell culture adaptations in virus stocks that can influence reagent generation and interpretation of a wide range of data including neutralization and drug efficacy. Overall, our work highlights the importance of this key motif in SARS-CoV-2 infection and pathogenesis. ARTICLE SUMMARY: A deletion of the furin cleavage site in SARS-CoV-2 amplifies replication in Vero cells, but attenuates replication in respiratory cells and pathogenesis in vivo. Loss of the furin site also reduces susceptibility to neutralization in vitro .

Movement of St. Louis encephalitis virus in the Western United States, 2014- 2018.

St. Louis encephalitis virus (SLEV) is a flavivirus that circulates in an enzootic cycle between birds and mosquitoes and can also infect humans to cause febrile disease and sometimes encephalitis. Although SLEV is endemic to the United States, no activity was detected in California during the years 2004 through 2014, despite continuous surveillance in mosquitoes and sentinel chickens. In 2015, SLEV-positive mosquito pools were detected in Maricopa County, Arizona, concurrent with an outbreak of human SLEV disease. SLEV-positive mosquito pools were also detected in southeastern California and Nevada in summer 2015. From 2016 to 2018, SLEV was detected in mosquito pools throughout southern and central California, Oregon, Idaho, and Texas. To understand genetic relatedness and geographic dispersal of SLEV in the western United States since 2015, we sequenced four historical genomes (3 from California and 1 from Louisiana) and 26 contemporary SLEV genomes from mosquito pools from locations across the western US. Bayesian phylogeographic approaches were then applied to map the recent spread of SLEV. Three routes of SLEV dispersal in the western United States were identified: Arizona to southern California, Arizona to Central California, and Arizona to all locations east of the Sierra Nevada mountains. Given the topography of the Western United States, these routes may have been limited by mountain ranges that influence the movement of avian reservoirs and mosquito vectors, which probably represents the primary mechanism of SLEV dispersal. Our analysis detected repeated SLEV introductions from Arizona into southern California and limited evidence of year-to-year persistence of genomes of the same ancestry. By contrast, genetic tracing suggests that all SLEV activity since 2015 in central California is the result of a single persistent SLEV introduction. The identification of natural barriers that influence SLEV dispersal enhances our understanding of arbovirus ecology in the western United States and may also support regional public health agencies in implementing more targeted vector mitigation efforts to protect their communities more effectively.

Measuring Alphavirus Fidelity Using Non-Infectious Virus Particles.

Mutations are incorporated into the genomes of RNA viruses at an optimal frequency and altering this precise frequency has been proposed as a strategy to create live-attenuated vaccines. However, determining the effect of specific mutations that alter fidelity has been difficult because of the rapid selection of the virus population during replication. By deleting residues of the structural polyprotein PE2 cleavage site, E3D56-59, in Venezuelan equine encephalitis virus (VEEV) TC-83 vaccine strain, non-infectious virus particles were used to assess the effect of single mutations on mutation frequency without the interference of selection that results from multiple replication cycles. Next-generation sequencing analysis revealed a significantly lower frequency of transversion mutations and overall mutation frequency for the fidelity mutants compared to VEEV TC-83 E3D56-59. We demonstrate that deletion of the PE2 cleavage site halts virus infection while making the virus particles available for downstream sequencing. The conservation of the site will allow the evaluation of suspected fidelity mutants across alphaviruses of medical importance.

Twenty years of West Nile virus spread and evolution in the Americas visualized by Nextstrain.

It has been 20 years since West Nile virus first emerged in the Americas, and since then, little progress has been made to control outbreaks caused by this virus. After its first detection in New York in 1999, West Nile virus quickly spread across the continent, causing an epidemic of human disease and massive bird die-offs. Now the virus has become endemic to the United States, where an estimated 7 million human infections have occurred, making it the leading mosquito-borne virus infection and the most common cause of viral encephalitis in the country. To bring new attention to one of the most important mosquito-borne viruses in the Americas, we provide an interactive review using Nextstrain: a visualization tool for real-time tracking of pathogen evolution (nextstrain.org/WNV/NA). Nextstrain utilizes a growing database of more than 2,000 West Nile virus genomes and harnesses the power of phylogenetics for students, educators, public health workers, and researchers to visualize key aspects of virus spread and evolution. Using Nextstrain, we use virus genomics to investigate the emergence of West Nile virus in the U S, followed by its rapid spread, evolution in a new environment, establishment of endemic transmission, and subsequent international spread. For each figure, we include a link to Nextstrain to allow the readers to directly interact with and explore the underlying data in new ways. We also provide a brief online narrative that parallels this review to further explain the data and highlight key epidemiological and evolutionary features (nextstrain.org/narratives/twenty-years-of-WNV). Mirroring the dynamic nature of outbreaks, the Nextstrain links provided within this paper are constantly updated as new West Nile virus genomes are shared publicly, helping to stay current with the research. Overall, our review showcases how genomics can track West Nile virus spread and evolution, as well as potentially uncover novel targeted control measures to help alleviate its public health burden.

Inter- and intra-lineage genetic diversity of wild-type Zika viruses reveals both common and distinctive nucleotide variants and clusters of genomic diversity.

Zika virus (ZIKV) strains belong to the East African, West African, and Asian/American phylogenetic lineages. RNA viruses, like ZIKV, exist as populations of genetically-related sequences whose heterogeneity may impact viral fitness, evolution, and virulence. Genetic diversity of representative ZIKVs from each lineage was examined using next generation sequencing (NGS) paired with downstream entropy and single nucleotide variant (SNV) analysis. Comparisons showed that inter-lineage diversity was statistically supported, while intra-lineage diversity. Intra-lineage diversity was significant for East but not West Africa strains. Furthermore, intra-lineage diversity for the Asian/American lineage was not supported for human serum isolates; however, a placenta isolate differed significantly. Relative in the pre-membrane/membrane (prM/M) gene of several ZIKV strains. Additionally, the East African lineage contained a greater number of synonymous SNVs, while a greater number of non-synonymous SNVs were identified for American strains. Further, inter-lineage SNVs were dispersed throughout the genome, whereas intra-lineage non-synonymous SNVs for Asian/American strains clustered within prM/M and NS1 gene. This comprehensive analysis of ZIKV genetic diversity provides a repository of SNV positions across lineages. We posit that increased non-synonymous SNV populations and increased relative genetic diversity of the prM/M and NS1 proteins provides more evidence for their role in ZIKV virulence and fitness.

Terrestrial Bird Migration and West Nile Virus Circulation, United States.

Host migration and emerging pathogens are strongly associated, especially with regard to zoonotic diseases. West Nile virus (WNV), a mosquitoborne pathogen capable of causing severe, sometimes fatal, neuroinvasive disease in humans, is maintained in highly mobile avian hosts. Using phylogeographic approaches, we investigated the relationship between WNV circulation in the United States and the flight paths of terrestrial birds. We demonstrated southward migration of WNV in the eastern flyway and northward migration in the central flyway, which is consistent with the looped flight paths of many terrestrial birds. We also identified 3 optimal locations for targeted WNV surveillance campaigns in the United States-Illinois, New York, and Texas. These results illustrate the value of multidisciplinary approaches to surveillance of infectious diseases, especially zoonotic diseases.

Genomic Characterization of Crimean-Congo Hemorrhagic Fever Virus in Hyalomma Tick from Spain, 2014.

Crimean-Congo hemorrhagic fever (CCHF) is a severe tick-borne disease caused by CCHF virus (CCHFV). Ticks in the genus Hyalomma are the main vectors and reservoirs of CCHFV. In Spain, CCHFV was first detected in Hyalomma ticks from Cáceres in 2010. Subsequently, two autochthonous CCHF cases were reported in August 2016. In this study, we describe the characterization of the CCHFV genome directly from Hyalomma lusitanicum collected in Cáceres in 2014. Phylogenetic analyses reveal a close relationship with clade III strains from West Africa, with an estimated divergence time of 50 years. The results of this work suggest that CCHFV has been circulating in Spain for some time, and most likely originated from West Africa.

Evidence for co-evolution of West Nile Virus and house sparrows in North America.

West Nile virus (WNV) has been maintained in North America in enzootic cycles between mosquitoes and birds since it was first described in North America in 1999. House sparrows (HOSPs; Passer domesticus) are a highly competent host for WNV that have contributed to the rapid spread of WNV across the U.S.; however, their competence has been evaluated primarily using an early WNV strain (NY99) that is no longer circulating. Herein, we report that the competence of wild HOSPs for the NY99 strain has decreased significantly over time, suggesting that HOSPs may have developed resistance to this early WNV strain. Moreover, recently isolated WNV strains generate higher peak viremias and mortality in contemporary HOSPs compared to NY99. These data indicate that opposing selective pressures in both the virus and avian host have resulted in a net increase in the level of host competence of North American HOSPs for currently circulating WNV strains.

Protective properties of vaccinia virus-based vaccines: skin scarification promotes a nonspecific immune response that protects against orthopoxvirus disease.

The process of vaccination introduced by Jenner generated immunity against smallpox and ultimately led to the eradication of the disease. Procedurally, in modern times, the virus is introduced into patients via a process called scarification, performed with a bifurcated needle containing a small amount of virus. What was unappreciated was the role that scarification itself plays in generating protective immunity. In rabbits, protection from lethal disease is induced by intradermal injection of vaccinia virus, whereas a protective response occurs within the first 2 min after scarification with or without virus, suggesting that the scarification process itself is a major contributor to immunoprotection. importance: These results show the importance of local nonspecific immunity in controlling poxvirus infections and indicate that the process of scarification should be critically considered during the development of vaccination protocols for other infectious agents.

Molecular epidemiology and evolution of West Nile virus in North America.

West Nile virus (WNV) was introduced to New York in 1999 and rapidly spread throughout North America and into parts of Central and South America. Displacement of the original New York (NY99) genotype by the North America/West Nile 2002 (NA/WN02) genotype occurred in 2002 with subsequent identification of a novel genotype in 2003 in isolates collected from the southwestern Unites States region (SW/WN03 genotype). Both genotypes co-circulate to date. Subsequent WNV surveillance studies have confirmed additional genotypes in the United States that have become extinct due to lack of a selective advantage or stochastic effect; however, the dynamic emergence, displacement, and extinction of multiple WNV genotypes in the US from 1999-2012 indicates the continued evolution of WNV in North America.

Continued evolution of West Nile virus, Houston, Texas, USA, 2002-2012.

We investigated the genetics and evolution of West Nile virus (WNV) since initial detection in the United States in 1999 on the basis of continual surveillance studies in the Houston, Texas, USA, metropolitan area (Harris County) as a surrogate model for WNV evolution on a national scale. Full-length genomic sequencing of 14 novel 2010-2012 WNV isolates collected from resident birds in Harris County demonstrates emergence of 4 independent genetic groups distinct from historical strains circulating in the greater Houston region since 2002. Phylogenetic and geospatial analyses of the 2012 WNV isolates indicate closer genetic relationship with 2003-2006 Harris County isolates than more recent 2007-2011 isolates. Inferred monophyletic relationships of these groups with several 2006-2009 northeastern US isolates supports potential introduction of a novel WNV strain in Texas since 2010. These results emphasize the need to maintain WNV surveillance activities to better understand WNV transmission dynamics in the United States.

The NYCBH vaccinia virus deleted for the innate immune evasion gene, E3L, protects rabbits against lethal challenge by rabbitpox virus.

Vaccinia virus deleted for the innate immune evasion gene, E3L, has been shown to be highly attenuated and yet induces a protective immune response against challenge by homologous virus in a mouse model. In this manuscript the NYCBH vaccinia virus vaccine strain was compared to NYCBH vaccinia virus deleted for E3L (NYCBHΔE3L) in a rabbitpox virus (RPV) challenge model. Upon scarification, both vaccines produced a desired skin lesion, although the lesion produced by NYCBHΔE3L was smaller. Both vaccines fully protected rabbits against lethal challenge by escalating doses of RPV, from 10LD(50) to 1000LD(50). A single dose of NYCBHΔE3L protected rabbits from weight loss, fever, and clinical symptoms following the lowest dose challenge of 10LD(50), however it allowed a moderate level of RPV replication at the challenge site, some spread to external skin and mucosal surfaces, and increased numbers of secondary lesions as compared to vaccination with NYCBH. Alternately, two doses of NYCBHΔE3L fully protected rabbits from weight loss, fever, and clinical symptoms, following challenge with 100-1000LD(50) RPV, and it prevented development of secondary lesions similar to protection seen with NYCBH. Finally, vaccination with either one or two doses of NYCBHΔE3L resulted in similar neutralizing antibody titers following RPV challenge as compared to titers obtained by vaccination with NYCBH. These results support the efficacy of the attenuated NYCBHΔE3L in protection against an orthologous poxvirus challenge.

Efficacy of CMX001 as a post exposure antiviral in New Zealand White rabbits infected with rabbitpox virus, a model for orthopoxvirus infections of humans.

CMX001, a lipophilic nucleotide analog formed by covalently linking 3-(hexdecyloxy)propan-1-ol to cidofovir (CDV), is being developed as a treatment for smallpox. In the absence of human cases of smallpox, new treatments must be tested for efficacy in animal models. Previously, we demonstrated the efficacy of CMX001 in protecting New Zealand White rabbits from mortality following intradermal infection with rabbitpox virus as a model for smallpox, monkeypox and for treatment of adverse reactions to smallpox vaccination. Here we extend these studies by exploring different dosing regimens and performing randomized, blinded, placebo-controlled studies. In addition, because rabbitpox virus can be transmitted via naturally generated aerosols (animal to animal transmission), we report on studies to test the efficacy of CMX001 in protecting rabbits from lethal rabbitpox virus disease when infection occurs by animal to animal transmission. In all cases, CMX001 treatment was initiated at the onset of observable lesions in the ears to model the use of CMX001 as a treatment for symptomatic smallpox. The results demonstrate that CMX001 is an effective treatment for symptomatic rabbitpox virus infection. The rabbitpox model has key similarities to human smallpox including an incubation period, generalized systemic disease, the occurrence of lesions which may be used as a trigger for initiating therapy, and natural animal to animal spread, making it an appropriate model.

Universal primers suitable to assess population dynamics reveal apparent mutually exclusive transcription of the Babesia bovis ves1alpha gene.

Babesia bovis is an intraerythrocytic hemoparasite of widespread distribution, which adversely affects livestock production in many regions of the world. This parasite establishes persistent infections of long duration, at least in part through rapid antigenic variation of the VESA1 protein on the infected-erythrocyte surface. To understand the dynamics of in vivo antigenic variation among the parasite population it is necessary to have sensitive and broadly applicable tools enabling monitoring of variation events in parasite antigen genes. To address this need for B. bovis, "universal" primers for the polymerase chain reaction have been designed for the ves1alpha gene, spanning from exon 2 to near the 3' end of cysteine-lysine-rich domain (CKRD) sequences in exon 3. These primers robustly amplified this segment, with minimal bias, from essentially the entire repertoire of full-length ves1alpha sequences in the B. bovis Mexico isolate genome, and are equivalently present in other isolates. On purified genomic DNA, this primer set can achieve a sensitivity of 10 genome equivalents or less. When applied to the amplification of cDNA derived from the B. bovis C9.1 clonal line evidence consistent with mutually exclusive transcription of the ves1alpha gene was obtained, concomitant with detection of numerous mutational events among members of the parasite population. These characteristics of the primers will facilitate the application of polymerase chain reaction-based methodologies to the study of B. bovis population and antigenic switching dynamics.